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ABSTRACT 

We discuss the phenomenology of SU(3) Skyrme and Nambu-Jona-Lasinio models. 
We show that while the Skyrme model, in its simplest version, is unable do reproduce 
the hyperon spectrum, the NJL model describes both hyperon and isospin splittings 
, with satisfactory accuracy. The difference between the two models lies in new 

■ anomalous moments of inertia, which vanish in the Skyrme model and get non- 

zero contribution from the valence quarks in the NJL model 
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1. Introduction 

Although QCD is now commonly accepted as the ultimate theory of strong in- 
teractions, the low energy properties of the hadronic states are usually calculated in 
terms of various effective models, which are believed to follow from QCD. Since there 
is however no rigorous direct derivation (see e.g.Ref.[l]), one has some freedom in 
choosing ones preferable model and preferable lagrangian. 
p i' In this talk we will discuss the phenomenology of the SU(3) Skyrme model^"^ and 
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the Nambu-Jona-Lasinio (NJL) model.^"^ In order to make the discussion transparent 
i) '. we will stick to the simplest versions of the models; the more refined versions are 

^ ' discussed here by other speakers. Our attention will be focused on the splittings both 

in SU(3) multiplets and in isospin multiplets. We will not consider absolute masses 
^ ■ which come out always too high. 

We will see that in the simplest version of the Skyrme model the pattern of split- 
tings is generically wrong: A— N= E — A = 2 (S — S) in contradiction with exper- 
imental situation where A— N^i S — S > S — A. To cure this disease one usually 
escapes to more complicated lagrangians and, at the same time, one apphes diago- 
nalization procedure of Yabu and Ando.^'' Although, as we will see on the example 
of the explicit perturbative calculation, this method shifts the spectrum in the right 
direction, one has to express a criticism against it, as it sums up an arbitrary subseries 
in the strange quark masses, neglecting other terms of the same order. 

In the semibosonised NJL model (see e.g. Ref. [11,12]), in which quarks interact via 
a self-consistent meson field, the spectrum is satisfactorily reproduced^'^ in the linear 
order in mg. In the NJL model the energy gets contribution from the valence and sea 
quarks. The classical part, i.e. the energy of the soliton, is exactly the same as in the 
two flavor case. The quantum corrections are calculated by adiabatical rotation of the 



*Talk presented at Siegen Workshop on Skyrem Model, October 1992, to be published in the 
proceedings 

'''Alexander von Humboldt Fellow, on leave of absence from the Institute of Physics, Jagellonian 
University, ul. Reymonta 4, 30-059 Krakow. Poland 



soliton resulting in a hamiltonian analogous to the one of the Skyrmion. A novelty is 
due to the mixed terms linear in the current quark mass and in the rotational velocity. 
These terms vanish in an ordinary Skyrme model; in the present model they get main 
contribution from the valence part. The resulting spectrum fits the data with a 10 % 
accuracy. At the same time the isospin splittings due to the rrid — mass difference 
are reproduced within experimental errors.^ Keeping in mind the simplicity of the 
starting lagrangian these results are surprisingly accurate. 

2. Gell-Mann— Okubo Mass Formulae 

Any dynamical model of light baryons has to reproduce and also explain the Gell- 
Mann-Okubo^^'^^ mass formulae which are derived assuming that the SU(3) breaking 
mass operator AM transforms like a F = 0, / = and /s = component of the octet 
tensor operator. Then, due to the Wigner-Eckhart theorem, matrix elements of this 
operator are given by: 

for the octet, and 

for the decuplet. The reduced matrix elements F, D and C are free constants. 
B = Y, I, Is for the baryon in question. The SU(3) Clebsch-Gordan coefficients can 
be written in terms of the diagonal SU(3) operators: 

AM|) = -^Y - ^(1 - + ^Y^), AM^' = -^Y. (3) 

The predictive power of Eq.(0) consists in the fact, that the number of free pa- 
rameters, which could in principle parametrize the mass splittings, is reduced from 3 
to 2 for the octet and to 1 for the decuplet. 

Equations (|) yield the relations: 

F = Ms- Mn, (4) 
-^D = J(Me - Ma) = ^(2Ms - Mh - Mn) = Mh + Mn - 2Ma, 

V5 2 6 

for the octet, and equal level spacing for the decuplet: 

-^C = Ms* - Ma = Afn. - Ms* = - M^*. (5) 
2v 2 

From these relations we can estimate the values of parameters F, D and C: 

F = 379, = 79 ± 17 and C = 415 ± 15 MeV. (6) 
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The resulting spectrum is presented in the first column of Fig.l. A small admixture 
of other operators like shifts the spectrum to the experimental position. 

The mass splittings of the baryons belonging to the same isospin multiplet consist 
of two parts: hadronic and electromagnetic/^ namely: 

ArriB = (AmB)h + {AmB)e- (7) 

If the isospin breaking is assumed to be driven by an octet isovector tensor operator 
corresponding to /a = then, in analogy to the Gell-Mann-Okubo mass formulae 
(Si), one gets: 

for the octet, and 

for the decuplet (normalization factors in front of the /, d and c are chosen for future 
convenience). Evaluating the SU(3) Clebsch-Gordan coefficients gives: 

{Am%\ = -i/l3 + rfYl3, {Amf\ = -^clg. (10) 

Electromagnetic part of the isospin splittings was estimated by Gasser and Leut- 
wyler^^ for the octet. Their estimate confirms a reasonable assumption that S~ — 
S"*" mass difference has no electromagnetic contribution. This assumption allows us 
to determine coefficient /, and also c for the decuplet, where no estimate of the 
electromagnetic part of S*~ — S*+ exists. Coefficient d can be determined from the 
hadronic part of the n-p mass difference.^ Altogether we get: 

/ = 12.11 ± 1.14 , ci= 1.73 ±0.38 and c = 6.6 ±1.0 MeV. (11) 

Symmetry considerations alone are not able to provide us with any relations be- 
tween the reduced matrix elements. Dynamical models, like Skyrme model or NJL 
model, make specific predictions for these constants. In the next sections we will 
calculate coefficients F, C, /, d, and c within the framework of the simplest 
version of the SU(3) Skyrme model and subsequently in the simplest version of the 
semibosonised SU(3) NJL model. 

3. SU(3) Skyrme Model 

Let us start by specifying the effective lagrangian proposed by Skyrme^^'^*^ and 
later generalized by Witten:^^'^° 

JdtLs^ = ^1 dtd^rTiid^U^d^U) + dtd'ri:i{[d^UU\d,UU^f) 

+ j dtd'ri:i{a{U + -2} + h{{U + U^)\s}] + N,T^f^i. (12) 
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Here is a number of colors, parameters and e, if taken from meson physics, 
are equal to ~ 186 MeV and e ~ 5.5 respectively. Fwz denotes the Wess-Zumino 
term^^'^^ and 

^K + m^), 6= ^K-m^). (13) 
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A time-independent Anzatz for the SU(3) matrix f/o takes the form of a hedgehog: 





cos P(r) + inr sin P(r) 















(14) 



where P(r) is a profile function. For the purpose of this talk we will take 



P(r) 



2 arctan 



(15) 



and calculate analytically all relevant quantities as functions of the soliton size tq. 
The energy of the solution (|1^) is then given by:^^ 
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with 



ml + 2ml 



K 



(412 MeV)^ 
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where Xq = eF^rQ. 

Minimizing Eq.(|l^) with respect to Xq we find in the chiral limit (/i^ = 0) Xq 



15/4 and M^i = 40.54 F^/e i.e. approximately 1370 MeV. This certainly unsatisfac- 
tory result (going off chiral limit shifts Md further up by about 280 MeV) is common 
for all chiral models which tend to overestimate the classical soliton mass. There is 
a hope that various other effects like gluonic corrections^^ or Casimir effect^^'^^ may 
bring this value down to less than 1 GeV. In what follows we will abandon the idea 
to fit the absolute masses, instead we will concentrate on the mass splittings. 
The splittings are calculated by rotating the static solution 



f/n 



A{t)UoA\t), 



where A G SU{3)/U{1), since [Ag, [/q] = 0. Therefore matrix A is defined up to 
a local U(l) factor h = exp(iA80),i.e. A and Ah are equivalent. This leads to a 
constraint which has to be imposed on the physical spectrum. Introducing 8 collective 
coordinates:^"^ 

AKt)j^A{t)='-j:X^n^; (19) 



a=l 
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one gets the following quantum mechanical lagrangian: 
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L = -M4P] + + + + Am[A, P], (20) 

^ i=l ^ k=4 

where Am [A, P] corresponds to the symmetry breaking piece. The moments of inertia 
can be again calculated analytically:^^ 

' 2^2 3 , ^ 107 



Ia = — rvr' 777 (63^0 + 25x0 



Ib = ^-^^(4x^ + 9x0) - (21) 

The standard quantization procedure^"^ leads to the following hamiltonian: 

Mb = Md + Hsvi2) + Hsvis) + i^br, (22) 

with 



C2(SU(2)r) C2(SU(3)l) - C2(SU(2)r) - g 

^SU(2) = , WSU(3) = ^ , (2dj 



where C2 denotes the quadratic Casimir operator of the right SU(2) symmetry corre- 
sponding to spin and of the left SU(3) corresponding to flavor. 

The wave function of the baryon state can be written as an SU(3) rotation matrix: 



MA) = v/dim(p, q) Dfl\A) = Jdun{p, q) (F, J, I, I D^^^'^\A) \ Fr, 5, -^3) , (24) 



where quantum numbers B and S denote now hypercharge, isospin and its third 
component, and right hypercharge, spin and its third component (with minus sign), 
respectively. SU(3) representations are labeled by (p, g), however not all p and q 
are allowed. The system is constrained since lagrangian (pO]) does not contain terms 
quadratic in the 8-th velocity. The constraint Yr, = A^c/3 selects the representations 
of triality zero: 

8, 10, To, 27, 35, 35, 64, . . . (25) 

for Nc = ?>■ The success of the model is the prediction that the lowest baryonic states 
belong to the octet and decuplet representations of SU(3). 
The symmetry breaking Hamiltonian: 

H^, = aD'il{A) (26) 

splits the hyperon spectrum (here representation (1,1) is denoted as (8) and index 8 
corresponds to F = and / = 0). Constant a is given as:^^ 

A/i2 3 7.7 GeV' 
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Fig. 1. Spectrum of light baryons. Dashed hnes represent experimental data, solid lines correspond 

to theoretical predictions described in the text. 
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where we have used 



A/i^ = -(m^ - ml) ^ (388 MeV)^ (28) 
3 

The first and probably the only good news about the breaking hamiltonian (^) 
is that it automatically fulfills the Gell-Mann-Okubo mass relation, since the left 
SU(3) Clebsch-Gordan coefficients coincide with the ones of Eqs([I|, H). The right C-G 
coefficients give specific predictions for the reduced matrix elements: 

F = --a, D = ^a, C = ^ a. (29) 

The spectrum however looks rather odd since the two ratios: 
F F 

— = = 2.24 (exp. 4.40), - = ^2 = 1.41 (exp. 0.91) (30) 

are far from their experimental values. If we choose a = —758 MeV in such a way that 
F is reproduced then we get spectrum shown in the second column of Fig.l. However 
a is not a free parameter. The 10-8 splitting (230 MeV) requires e^F^^ = 16.41 GeV. 
Through Eqs(|T], |3) we get that 1/Ib = 405 MeV and a ^ -470 MeV, approximately 
40 % smaller than the value required to fit F. One could in principle enlarge a by 
tuning the kaon mass, but the problem of Eq. (|30D remains. 

Guadagnini^ in the first paper on the SU(3) quantization of the Skyrme model 
proposed to add by hand a term proportional to the hyphercharge: 

HG = aD^^l{A)+pY. (31) 

(3 is a free parameter and the reduced matrix elements read: 

F = -^a-2p, D = -^a, C = -^(a + 8/3). (32) 

An excellent fit to the data (column 3 in Fig.l) is obtained with a = —382 MeV and 
(3 = —94 MeV. This time a is smaller than the value required by 10-8 splitting, but 
only by about 20 %. 

One way to generate effective (3 is to calculate the O(a^) correction to baryon 
energy (this approach for not too large a is in fact equivalent to the diagonalization 
procedure of Yabu and Ando^'' which sums up the perturbative series in a): 
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Constants F, D and C can be immediately read off from Eqs(^), moreover new 
terms not present in the original Gell-Mann-Okubo mass formula (|^) are generated. 
Numerically they are small and we will neglect them in the present discussion. 

In order to fit F and D one needs a ~ —677 MeV and 1/Ib = 378 MeV; with 
these values one gets a bit too small C = 327 MeV. These values of a and Ib are 
again larger than the ones required by 10-8 splitting. 

This was the status of the SU(3) Skyrme model in the mid eighties. A philosophical 
question whether one should trust perturbative expansion in the strange quark mass 
arose. Some people decided to work within an approach which breaks the SU(3) 
flavor symmetry already at the level of the Anzatz for the Uq field, ^^'^^ others tried to 
push up a and Ib either by enlarging the original Skyrme lagrangian or by employing 
another model (like NJL model discussed in the next section). 

4. Semibosonised SU(3) Nambu— Jona-Lasinio Model 

It was already shown in the talk of Alkofer^^ that the solitonic solutions of the 
semibosonised NJL model are studied in terms of an effective action: 



S^ff = -Splog(i^ - m - MU^'). 



(34) 



This time, however, U is an SU(3) matrix of Eq.(|T^. M is the constituent quark 
mass, which is in fact the only free parameter of the model. The bare quark mass 
matrix can be written in a form: 



m = /io Ao - /Us As - /is A3, 
where Aj are Gell-Mann SU(3) matrices (Aq = v2/3 1) and 



(35) 



1^0 



1 



1 



/is 



(2ms-mu-md), /is = ("^d - "^u)- 



(36) 

The energy of the soliton consists of two parts: the energy of the continuum, i.e. 
the energy corresponding to the effective action (|3^) , and the energy of the valence 
level. In what follows, for simplicity, we confine our discussion to the continuum part 
only, but one always has to remember that the pertinent valence contribution has to 
be added. 

The effective action (Q) can be rewritten in terms of the Euclidean spectral 
represent at ion : ^ ' ^ 



eflf 



-NrT 



duj 
2^ 



Tr log {ioj + H) 



1 + 



1 



iuj + H 



(37) 



where H is the hermitean static hamiltonian: H = 74(7i5j + MUq), and appropriate 
regularization is understood. The static soliton solution for H reduces to the one 
found in the SU(2) case. Formula (^) is already written in a form ready to be 
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expanded in a power series in m and in generalized velocities Q. Let us for the 
moment forget about the mass matrix m and expand (0) in powers of Q. We get 
(back in Minkowski metric) familiar lagrangian (^) with: 



tab 



2ti 



-Tr 



1 



iu! + H iuj + H 



A. 



lA^ab for a, 6 = 1...3 
IsSab for a, 6 = 4. ..7. 
for a, 6 = 8 



(38) 



The above expression is assumed to be properly regularized and the full moments of 
inertia have also a valence part. 

The quantization proceeds exactly as in the case of the Skyrme model, and as 
a result one arrives at the hamiltonian ([2^). 

The novelty comes from the expansion in powers of the rotated matrix m: 



l(8) 



(8) 



l(8)> 



(39) 



where constant a 



1 



iu + H 



74 Aa 



-\/6 a for a = 
a for a = 8 
for a = 1...7 



(40) 



is related to the pion-nucleon sigma term S = 3/2{m^ + rridjcr and anomalous tensor 
Kab is defined as: 



Kab = ^ — / — Tr 
4 J 2n 



1 



iuo + H 



74AaT 



iuj + H 



KA^ab for a,b = 1...3 
KB5ab for a, 6 = 4... 7 . (41) 
for a, 6 = 8 



We call Kab anomalous since it comes from the imaginary part of the effective action, 
which is related to anomaly, and as such does not require regularization. In fact Kab 
gets contribution almost entirely from the valence level. 

The quantized hamiltonian i^br is not as simple as in the Skyrme model and reads: 



where 



br 



Dg\A) + /38 Y + 78 T.D'sJ{A)S 



(8), 



a=l 



a=l 



-a + 



K, 



Pi = -VS^ Hi, 7j 



Ka K 



B 



IB 



Index i = 8 corresponds to F = 0, / = and index i = 3 to F = 0, / = 1, Is 



(42) 



(43) 



0. 
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We see that ( ^2|) contains a term proportional to Y which naturally arises in this 
model, moreover there is another term proportional to the product of a D function 
and the spin operator. 

We adopt the following numerical procedure: first we find the solitonic solution 
for a range of constituent masses M, then we find the optimal value of M which 
reproduces the 10-8 splitting due to the rotational hamiltonian -ffsu(2) (p3D. It turns 
out®' that M = 390 MeV and the corresponding moments of inertia take the following 
values: 1/Ia = 157, 1/Ib = 234, 1/Ka = 469 and 1/Kb = 704 MeV, and a = 
3.07.^^ We will see that not only hyperon splittings but also isospin splittings are well 
reproduced. To this end let us let us define the following quantities: 

^ = a + 2— + —, 7 = a + 2— + 5— , 6 = a + 2—-3—. (44) 

Kb Ka Kb Ka Kb Ka 

Then we get: 

F = -09 777,8, D = — ^5ms, C = — -='yms, (45) 

2^ " 2^5 ' 2^2 ' ^ ' 

3 3 3 

/ = - V2 (r?7d - mu), (i = — 5(md -mu), c = - 7 (ma - mu). (46) 

4 /U o 

In Figs 2 and 3 we plot the splitting constants of Eqs(^, ^ ) as functions of rus 
and md — respectively, together with the error bars corresponding to Eqs(P, |Tl|). 
For r?7s ~ 175 MeV the N-S splitting (i.e. constant F) is reproduced. For the same 
value of rris constant D corresponding to E-A splitting is overestimated by 35 MeV, 
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MeVl 



Fig. 3. Constants /, d and C. Solid lines represent model predictions as functions of md — rriu 
mass difference. Dashed line correspond to the error bars of Eq.(ll). 



whereas constant C is underestimated by 60 MeV. The resulting spectrum is shown 
in column 4 of Fig.l. On the other hand, the isospin breaking constants /, d and c 
are reproduced within experimental errors for itla — ^ 3.5 MeV (see Fig. 3). 

The present model makes specific prediction for the ratio of hadronic to isospin 
breaking constants: 

f d c 

4; (2.99 ±0.3) = VS— (4.90 ±2.15) = 72- (2.25 ±0.42) (xlO^^), (47) 

where the numbers in brackets correspond to the experimental values of Eqs(P, [Tl|). 
Certainly the central values are fairly scattered. We would like to offer the following 
explanation of this discrepancy. The isospin splittings are proportional to a tiny 
parameter, namely ma — mu, and therefore the first order of the perturbation theory 
is legitimate. On the contrary, for the hyperon splittings which are proportional to the 
much larger parameter, namely strange quark mass, one may expect some corrections 
from the higher order of the perturbative expansion in rris. Indeed, already the second 
order brings the splittings to their experimental values with an accuracy of a few 
MeV. Last column in Fig.l shows the spectrum obtained by applying the Yabu Ando 
procedure to the present model. 

To summarize: we have studied the symmetry breaking effects due to the quark 
masses in the Skyrme model and in the solitonic sector of the semibosonised NJL 
model. In the simplest version of the Skyrme model the splitting operator is too 
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simple to account for hyperon splittings. On the contrary, in the simplest version 
of the NJL model, a satisfactory description of the hadronic mass spectrum includ- 
ing both the hyperon and the isospin splittings was found. The new terms in the 
splitting operator, which are not present in the Skyrme model, are proportional to 
the anomalous moments of inertia which get the main contribution from the valence 
level. The absolute masses are too big, but there exist several mechanisms which 
may bring them down, namely gluon corrections,^^ rotational and translational band 
subtraction^ and Casimir energies of quantum fluctuations. ^^'^^ 
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